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I. INTRODUCTION 



Multiple parton interactions are one of the most com- 
mon, yet poorly understood [l|, phenomena at the LHC. 
It is usually assumed that only one hard interaction takes 
place per hadron-hadron collision together with multiple 
soft interactions. This assumption is typically justified 
on the grounds that the probability of a hard parton- 
parton interaction in a collision is very small. Thus the 
probability of having two or more hard interactions in a 
collision is strongly suppressed in comparison with the 
single interaction probability. Nevertheless hadron col- 
lisions in which two (or more) distinct pairs of partons 
hard scatter are possible. 

The presence of such multiple parton interactions in 
high energy hadron-hadron collisions has been convinc- 
ingly demonstrated by the AFS @, UA2 @ and CDF Q 
collaborations using events with the four-jet final state 
and later by the CDF Collaboration Q and quite re- 
cently by the DO Collaboration Q using events with the 
7 + 3 jets final state. The possibility of observing two 
separate hard collisions was proposed long ago and their 
theoretical study goes back to the early days of the par- 
ton model LZripl with subsequent extension to perturba- 
tive QCD [l(H23. 

A greater rate of events containing multiple hard in- 
teractions is anticipated at the LHC with respect to the 
experiments mentioned above due to the much higher lu- 
minosity and greater energy of the LHC. Moreover the 
products from multi ple interactions will represent an im- 
portant background [29r{32| to signals from the Higgs and 
other interesting processes and certain types of multiple 
interactions will have distinctive signature j33T - [3r| facili- 
tating a detailed investigation of these processes experi- 
mentally. Therefore it is extremely important to combine 
theoretical efforts in order to achieve a better description 
of multiple interactions, in particular, double scattering 
which will be the dominant multiple scattering mode at 
the LHC. 

With the only assumption of factorization of the two 
hard parton processes A and B the inclusive cross section 
of a double parton scattering process in a hadron collision 



may be written as [Io|, [H, IH[ 




xof l {x2-, x 2 )T M (.x 1 , x 2 ,b; Q 2 , Q 2 2 )dx\dx 2 dx x dx 2 d 2 b. (1) 

Here Fy (xi,x 2 ,b;Ql,Q 2 ) are the double parton distribu- 
tion functions, depending on the longitudinal momentum 
fractions x\ and x 2 and on the relative transverse dis- 
tance b of the two parton undergoing the hard processes 
A and B at the scales Q\ and Q 2 . a^ k and are the 
parton- level subprocess cross sections. The factor m/2 
is a consequence of the symmetry of the expression for 
interchanging parton species i and j. m = 1 if A = B 
and m = 2 otherwise. 

It is often assumed that the double parton distribution 
functions may be decomposed in terms of longitudinal 
and transverse components as follows: 

r« (a*, x 2l b; Ql Ql) = D%( Xl ,x 2 ;Ql Q 2 2 )F l3 (b). (2) 

Fij{b) represents the parton pair density in transverse 
space and is not calculated in the framework of pertur- 
bative QCD. Existing models typically use Gaussian or 
exponential form to describe the Fij (b) , or a sum of Gaus- 
sian/exponential terms [22I I37I. 

D 1 ? [x\, x 2 ; Qi, Q 2 ) is interpreted as the inclusive prob- 
ability of finding a parton i with the longitudinal mo- 
mentum traction x\ at scale Qi and a parton j with 
the longitudinal momentum traction x 2 at scale Q 2 in a 
hadron h. The behavior of these distribution functions 
D l ^(xi,x 2 ;Q 2 ) = D^(xi,x 2 ;Q 2 ,Q 2 ) with the two hard 
scale set equal were investigated in Refs. [3^, H^] in the 
leading logarithm approximation of perturbative QCD. 
The equations dictating the scaling violations (i.e. Q 2 
dependence) of the double parton distributions were de- 
rived therein. A crucial prediction {40TJ42j of this is that 
even if the double parton distribution functions may be 
taken to be equal to products of single distribution func- 
tions at one scale, then at any other scale the double 
distribution functions will deviate from factorized forms 
used often in the current analysis. Quite recently the 
equations of Refs. [H, HU have been numerically inte- 
grated [25[ to produce a set of publicity available double 
parton distribution function grids covering the ranges 
1(T 6 < an < 1, 10" 6 < x 2 < 1, and 1 < Q 2 < 10 9 
GeV 2 . Thus the correlations induced by QCD evolution 
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may be evaluated already under condition of real exper- 
iments and their possible manifestation was discussed in 
Refs. [lil, US E3 • In this case the behavior of double par- 
ton distribution functions is extremely important in the 
kinematical range of relatively small longitudinal momen- 
tum fractions which is difficult to investigate numerically 
with a good accuracy and predictive enough power. 

The main purpose of the present paper is to study the 
asymptotic behavior of double parton distribution func- 
tions in the vicinity of the kinematical boundaries analyt- 
ically. The paper is organized as follows. Section|TI]is de- 
voted to what may be extracted from perturbative QCD 
theory on the asymptotic behavior of the two-parton dis- 
tribution functions at the parton level. The possible phe- 
nomenological issues for parton distributions in a hadron 
are discussed in Sec. IIIII We summarize and conclude in 
Sec. El 



DOUBLE PARTON DISTRIBUTIONS IN 
THE LEADING LOGARITHM 
APPROXIMATION AT PARTON LEVEL 



II. 



The analysis [45M47| of hard processes, ep scattering 
and e + e~ annihilation, in the leading logarithm approx- 
imation of perturbative QCD admits a simple interpre- 
tation [46( | in the framework of the parton model with 
a variable cutoff parameter A ~ Q with respect to 
the transverse momenta. The dependence of the multi- 
parton distribution and fragmentation functions on the 
value of this cutoff parameter is determined by the evo- 
lution equations. Lipatov suggested a general elegant 
method |46j of obtaining these equations in any renor- 
malizablc quantum field theory. As an evolution vari- 
able one uses the value of hard scale (the transfer mo- 
mentum squared Q 2 , most commonly) or its logarithm 
£ = ln(Q 2 //i 2 ) or double logarithm 



t = 



1 

2^(3 



In 



-V) 




1 

2^5 



In 



) 



QCD 



ln(- 



-) 



QCD J 

(3) 

where (3 = (33 — 2n^)/127r in QCD, g(ii 2 ) is the running 
coupling constant at the reference scale [i 2 , n/ is the 
number of active flavors, Aqcd is the dimensional QCD 
parameter. 

With the evolution variable t (taking into account the 
behavior of the running coupling constant in one l oop 
approximation explicitly) the DGLAP equations [45l448| 
for single distributions read more easily 



dDj(x,t) 
dt 



E 



^ D ((x',t)P r ^ 



(4) 



They describe the scaling violation of the parton distribu- 
tions Dj [x, t) inside a dressed quark or gluon (i, j = q/g). 
It is interesting that expression for the kernels P in Lipa- 
tov method already includes a rcgularization at x — > x' , 



which was introduced in Ref. [48[ afterwards based on 
the momentum conservation. 

This method allows one to obtain also the gener- 
alized DGLAP equations for two-parton distributions 
Df 1J2 (x\, X2, t), representing the probability that in a 
dressed constituent i one finds two bare partons of types 
ji and J2 with the given momentum fractions x\ and x%, 
namely (see Refs. (H, [H[ for details): 



dD^ 2 (x 1 ,x 2 ,t) 



dt 

i—x 2 



(5) 



E 

ii' 



Xl 



jDf ^(x 1 ',x 2 ,t)P jl ^ jl 



xi' 



1 — CCl 



E J ^WOJWte) 

JV X2 \ J 



where the kernel 



-Pr 



Xl 



X! +x 2 ^ 3l32y x 1 +x 2 ' 



(0) 



is defined without <5-function rcgularization. 

It is readily verified by direct substitution that the so- 
lution of Eq. ([S]) can be written via the convolution of 
single distributions [3^, Ht| 

D^ 2 (x u x 2 ,t) (7) 

t 1 — X2 1 — zi 




-2 



This coincides with the jet calculus rules [49( proposed 
originally for the fragmentation functions and is the gen- 
eralization of the well-known Gribov-Lipatov relation in- 
stalled for single functions [45l - [47j (the distribution of 
bare partons inside a dressed constituent is identical to 
the distribution of dressed constituents in the fragmen- 
tation of a bare parton in the leading logarithm approx- 
imation). Nevertheless in a phenomenological employ- 
ment the direct numerical integration of Eq. ([S]) occurs 
more effective [25| than the explicit solution (JT]) because 
of the singular i5-like initial conditions of single distribu- 
tions Dj (x, t) (the Green's functions) incoming in this 
solution. 

Eqs. (f5]) and ([5]) are explicitly solved by introducing 
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the Mcllin transforms 



M/'(n,t) = / dxx n D j h (x,t), 
o 



(8) 



= Jd Xl dx 2 0(l - Xl - x 2 )x n 1 1 X , 2 l2 D j h lj2 (x 1 ,X2,t) (9) 



which lead to a system of ordinary linear-differential 
equations at the first order: 

dMf{n,t)/dt = Y,M{ {n,t)P r ^{n), (10) 



dM 3 > j2 (n u n 2 ,t)/dt 

= J2Mt' j2 (ni,n 2 ,t)P h ^ h ( ni ) 
h' 

+ J2Mt h '(n 1 ,n 2 ,t)P j2 ^ j2 (n 2 ) 

37.' 

+ ^2 M i ( n l + n 2,t)Pj'^j!j2( n l, n 2), (11) 



where 

l 

P^ n (n) = J x n P i ^ il (x)dx, (12) 
o 

l 

Pi^ ilh {n l7 n 2 ) = J x ni (1 - x) n2 P t ^ lll2 (x)dx. (13) 
b 

In order to obtain the distributions in x representation 
an inverse Mellin transformation must be performed 

xD{(x,t)= J ^x-»M/'(n,t), (14) 

where the integration runs along the imaginary axis to 
the right from all n singularities. This can be done in a 
general case only numerically. However the asymptotic 
behavior can be estimated in some interesting and partic- 
ularly simple limits in the technique under consideration. 



A. Spectra at finite x 

At first let us look at single parton distribution func- 
tions near x = 1. These can be obtained through the 
study of large moments (which clearly emp hasize x ~ 1). 
One obtains behaviors of the type [43, |49| : 



Dl{x,t)\ x ^~Cl{x,t){l-x) 



2 C F t+bi 



(15) 



where C F = (7V C 2 - 1)/2A C = 4/3, N c = 3 and Cj(x,t) 
depend logarithmically upon (1 — x) and in a compli- 
cated but known way upon t, whereas bij are just some 



integer numbers. The values of Cj (x, £), bij are given in 
Table 4 of Ref. [Hj](NP) for all cases. For illustration we 
write here the asymptotic expression for the distribution 
of valence quarks (the probability of finding in a quark a 
quark of the same species) at 1 — x <C 1 H3, S^] : 



D val (x, t) = cxp ((3/2 - 2 lE )C F t) 



(1-x) 



2 C F t-l 



T{2C F t) 



(16) 



where ^ F = 0.5772... is the Euler's constant. Indeed, at 
large n the kernel 



P q ^ q (n) ~ (3/2 - 2 lE - 21nn)C F , 



(17) 



and an inverse Mellin transformation (|T4|) can be per- 
formed analytically using the exact solution 



M val (n,i) = exp(P,_>,(n)t) 



(18) 



in n representation. The characteristic moment being 
significant in integral (|14[) is defined by the saddle point 
equation 



no = 2C F t/{\ - x), n > 1. 



(19) 



For the two-parton distribution functions the limit 
x\ + x 2 — > 1 can be investigated by looking into large 
moments: 



1 



dxxdx^X! +x 2 ) n x 1 x 2 Dl 1 n (x u x 2 ,t). (20) 







The result is a behavior of the type [49( : 

Df h (n, i 2) t) ~ A? h (1 - Xl - x 2 f CFt+Si ^ , (21) 

where j4| U2 have at most a logarithmic dependence on 
(1 — xi — x 2 ). The exponents 6 j. , are some integer 
numbers given in Table 5 of Ref. |49j(NP). 
The analog of a doublc-Rcggc limit 



l-xi<l, 1 — xi — x 2 <C 1 — xi 



(22) 



can be done by studying double moments Mf 132 (m, n 2l t) 
for ni 3> n 2 3> 1. By inverting the Mellin transform one 
then finds 



D^( Xl ,x 2 ,t) 

~ H(t)(l - xi) ki ^ (1 - xi - x 2 ) 2Cpt+h 



up to logarithmic terms. The exponents k 



1,3132 



(23) 



and 



hi.ii j ? can be computed and are given in Table 6 of 
Ref. gg|(NP). 

The obtained asymptotic behaviors ([TBI |2"T1 |2"3"|) admit 
an amusing analogy with Regge theory as it was pointed 
out in Ref. [49} . This helps one also to understand better 
why the phase space factor in the improved ansatz [25[ 
for the initial conditions of parton distributions at the 
hadron level must include the factors (1 — x\ — x 2 ), (1 — 
xi),(l — x 2 ) with the exponents depending on parton 
types in order to satisfy the momentum and number sum 
rules. 
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B. Spectra at small x 

The increase in the average number of decays of each 
parton with increasing t leads, as can be readily under- 
stood, to a large increase in the number of fields, glu- 
ons and sea quarks, with small longitudinal momentum 
fractions. Technically, integral (| 14[) when x <C 1 is de- 
termined by the first singularity with respect to the n of 
the moments M\ (n, t), which is related to the pole of the 
larger eigenvalue A+(n) of the system (fTCTj) when n = 0: 



n 



n.< 1, 



(24) 



where a = ^-N c + \n f (l - 2C F /N C ). By inverting the 
Mellin transform with the approximation (|24"|) one then 
obtains the probability of finding a gluon in a gluon: 

xD 9 g (x,t) = / — exp[nln(l/a;) + 2N c t/n- at] 
= 4,N c texp[-at]Ii(v)/v 

~4:N c tv- 3/2 exp[v-at]/V2^, (25) 



where v = y^8N c t In (1/x) and 1\ is the usual modi- 
fied Bessel function. The last approximate equality in 
Eq. (f2"5")) is obtained by the saddle point method with the 
saddle point equation 



n = y / 2N c t/\n(l/x), n <C 1. 



(26) 



The probabilities of finding any other partons in a gluon 
or a quark are evaluated in the regime of small x by a 
similar way and can be found in Refs. j47l.l49j. 

Due to the generalized Gribov-Lipatov relation men- 
tioned above in the limit of one slow (x± ~ 0) and one 
fast (x 2 = finite) parton the two-parton distribution func- 
tions become uncorrelated the same as the two-parton 
fragmentation functions for which this property was es- 
tablished in Ref. [4!|. The factorization property is sig- 
nificant in a phenomenological employment and its proof 
is worth to be reviewed once more because of a very old 
affair. 

The exact solution for single distributions in the mo- 
ments representation can be written symbolically in a 
matrix form 



M{{n,t) = (cxpP(n)^' 



(27) 



The solution for the two-parton distribution functions in 
this representation reads 



M^ 2 (n 1} n 2} t) 



Z 

= J2 /^'(expPK+n 2 )i')^[P^ VjV ( 



ni,n 2 )} 



Going further through a Laplace transform to the vari- 
able E one finds 



Mf j2 (m,n 2 ,E) = / dte- m Mi lj2 (ni,n 2 ,t) 



Yl (p — UTzr i ~ )j[ J j->jyjv("i»"a)] 



x( 



f -f. ,E-P{m + m) 

1 yih 

E-P( ni ) -P(n 2 ) h ' j2 ' 



(29) 



Let us consider now the correlation function and its mo- 
ments: 

Cf^(x 1 ,x 2l t) 

= D^ 2 (x x ,x 2 ,t) - D{ 1 (x 1 ,t)Di 2 (x 2 ,t), (30) 



Cf^{e^E) 



(31) 



x( 



E - P(e) - P(n) 



\hh _ ( \hh 

hi' h' K E-P(e)-P(n) hl ' 



x(expP(n 1 )(t - 0)5i/(expP(n a )(< - ?))%,. (28) 



where e = + in the moment relative to the slow parton 
and n = 0(1). In Eq. (|3"Tj) various factors are divergent 
since both P(e) and Pj^j 1 'j 2 > (e, n) contain divergent ma- 
trix elements. Since, however, 

Pwh'k, n) ~ P jfi (e) P(e + n) ~ P(n) + 0(e) (32) 

one finds that Cf l3Q '(e,n,E) is NOT divergent. In other 
words the divergent part in Df lj2 (xi, x 2 ,t) at x\ — > 
and finite x 2 factorizes into D J i 1 (x x ,t) times D 02 (x 2 ,i). 
This result can be generalized to m-parton distribution 
functions, i.e., when one of the partons becomes wee its 
divergence factors out. 



III. DOUBLE PARTON DISTRIBUTIONS IN 
HADRONS 

Of course, it is interesting to find out the phenomeno- 
logical issue of this parton level consideration. This can 
be done within the well-known factorization of soft and 
hard stages (physics of short and long distances) ■ As 
a result the equations ((H [TUl QJ} describe the evolution 
of parton distributions in a hadron with t (Q 2 ), if one re- 
places the index i by index h only. However, the initial 
conditions for new equations at t = (Q 2 = /j 2 ) are 
unknown a priori and must be introduced phenomeno- 
logically or must be extracted from experiments or some 
models dealing with physics of long distances [at the par- 
ton level: D{(x,t = 0) = 6 tJ S(x - 1); Df J2 (x u x 2 ,t = 
0) = 0]. Nevertheless the solutions of the generalized 
DGLAP evolution equations with the given initial con- 
ditions may be written as before via the convolution of 
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single distributions and in the moments representation 
they read 

Mi lh (n lt n 2 ,t) = M*g*(n lt n2,t) 
+ E M^' A '(n 1 ,n 2) 0)A^,(ni > t)M?,(n 2> t), (33) 



h'h' 



where 



Y,M t h (n 1 +n 2 ,0)Mi ll2 (n 1 ,n 2 ,t) 



(34) 



are determined by the first singularity with respect to the 
ri\ of the moments M d ^ n (ni, n 2 , t) with n 2 = 0(1). In 
this case in accordance with the consideration above the 
perturbatively calculated particular solutions ([3"5]) of the 
parton level factorizes into M^^ni^t) times M\ 2 (n 2 ,t). 
Thereafter the general solutions (|53"|) read 

M£ h ( ni ,n 2 ,t) (37) 
= J2 M i(ni + n 2 , 0)Mf (m,t)Mi 2 (n 2 ,t) 

i 

+ E Mi 1 \n 1 ,0)Ml 2 \n 2 ,0)Mjl,(n u t)M^(n 2 ,t), 



are the particular solutions of the complete equations 
with the zero initial conditions at the hadron level and 

Mi l32 ( ni ,n 2 ,t) = J2 I ' dt'Mi{n x + n 2 ,t') 

Pw^inun^Mfrinut- t')M%(n 2 ,t - t') (35) 

are the particular solutions of the complete equations 
with the zero initial conditions at the parton level. While 
the second term in the expression (|33[) represents the so- 
lutions of the homogeneous evolution equations with the 
given initial conditions Mj^ 32 (ni,H2,0). These nonper- 
turbativc unknown two-parton initial conditions are just 
the reckoning for the unsolved confinement problem. If 
one assumes that in the moments representation there is 
the approximate factorization of these initial conditions 

M^'tm.na.O) ~ M£'( ni ,0)Mi 2 '(n 2 ,0), (36) 

then the solutions of the homogeneous evolution equa- 
tions will be factorized in x representation. 

The exact solutions (|33p in the moments representation 
together with some additional reasonable assumptions on 
the initial conditions allow one to extract the asymptotic 
behavior of the distribution functions already in hadrons 
near the kincmatical boundaries as it has be done for the 
functions of the parton level. Since the kinematical range 
of relatively small longitudinal momentum fractions is 
extremely important under condition of real experiments 
we focus on the asymptotic behavior of two-parton dis- 
tribution functions just in this range. 

However some interesting results at finite x can be 
found in Ref. [HT| f° r two-particle fragmentation func- 
tions and can be extended to two-parton distribution 
functions taking into account the difference between the 
solutions of homogeneous evolution equations for the 
fragmentation and distribution functions (52j. Thus the 
asymptotic behaviors like (f^Tl can be obtained at 
the hadron level also assuming that the initial conditions 
can be evaluated in the saddle points which were defined 
at the parton level or / and fixing the asymptotic form of 
initial conditions near the kinematical boundary. 

In the important limit of one slow (x\ ~ 0) and one fast 
(x 2 = finite) parton the two-parton distribution functions 



the natural approximate factorization of the initial con- 
ditions in the form (|36[) being supposed. 
Since in the limit e — > + [43, |49| 



Mf(e,t) 
M«(e,t) 



No 
C F 
Nc 



M|(e,t)=exp(A + (e)t), 



(38) 



6JV„ 



exp(A+(e)i),(39) 



then the ratio of the particular solutions (correlations 
induced by evolution) to the solutions of homogeneous 
equations (giving approximate factorization) reduces to 



M 9 h £je,n 2 ,t) 



(40) 



J2Mi(e + n 2 ,0)Mf(e,t)Mt(n 2 ,t) 

_ i 

= E W (*> WL> (*> *) E K 2 ' M<K *) 

h ' h ' 

M 9 h (n 2 , 0)M*(n 2> t) + % £ M«(n 2 , 0)Mf (n 2 ,t) 

_ 9 

M(e,0) + ^E^(e,0)]M*(n a ,t) 
<? 

and the universal divergent factor exp (A + (e)i) is can- 
celled out. One should note that this cancellation takes 
place also in the general case without the supposed fac- 
torization of initial conditions in the form (|36[) . Therefore 
the ratio (|4T)|) admits the finite limit at e — > + provided 
that the initial number of quarks and gluons in a hadron 
is finite, i.e. M^(0, 0) and M^(0, 0) are not singular and 
can be evaluated in the saddle point (|26|) . Without a fac- 
tor C f /N c before the quark sum the numerator in the ra- 
tio under consideration (|40[) is simply equal to M^ 2 (n 2 , t) . 
Therefore the suppression of correlations induced by evo- 
lution with respect to the factorization component will be 
mainly determined by the gluon and quark multiplicities 
in the initial conditions (with a possible correction re- 
lated to the estimation of zero moments in the the saddle 
point): 



D±Jx u x 2 ,t) 



xi->0 



M«(0,0) + ^£A^(0,0)' 



(41) 
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The similar result is obtained if one slow quark is con- 
sidered instead of a slow gluon. 

Thus in the important limit of one slow [x\ ~ 0) and 
one fast (x2 = finite) parton the two-parton distribution 
functions possess practically the factorization property. 
The additional contribution induced by evolution being 
suppressed by the initial gluon and quark multiplicities in 
comparison with the "genuine" factorization component 
(the solution of homogeneous equation). 

IV. CONCLUSIONS 

The asymptotic behavior of the two-parton distribu- 
tion functions has been investigated near the kinematical 
boundaries. At x close to 1 these functions include the 
factors (1 — xi — X2), (1 — (1 — x 2 ) with the exponents 
depending on parton types. These exponents are known 
at the parton level and can be calculated in principle at 
the hadron level fixing the asymptotic form of initial con- 
ditions near this kinematical boundary. The importance 
of the factors under consideration was found in Ref. [25| 
at the numerical integration of evolution equations in x 
representation taking into account the momentum and 
number sum rules. 

The two-parton distribution functions become prac- 
tically uncorrelatcd in the kinematical range of rela- 
tively small longitudinal momentum fractions which is 
extremely important under condition of real experiments. 
The additional "factorization" contribution induced by 



evolution being suppressed by the initial gluon and quark 
multiplicities in comparison with the "genuine" factoriza- 
tion component (the solution of homogeneous equation) 
in the case of one slow {x\ ~ 0) and one fast (X2 = finite) 
parton. This revealed property explains why the indica- 
tions from the experimental observation of double parton 
scattering at CDF Q are not in favor of correlation ef- 
fects in longitudinal momentum fractions. It is extremely 
difficult to select events experimentally with double par- 
ton scattering in which both partons are enough fast and 
therefore are strongly correlated in longitudinal momen- 
tum fractions. In this connection it is interesting to study 
the double parton distribution functions beyond the lead- 
ing logarithm approximation over Q 2 (for instance in the 
BFKL regime |53l - l56l | Q 2 = const and \n(l/x) — > 00 
with the possible nontrivial nonfactorization behavior of 
the two-parton distribution functions just at small x). 
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